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ABSTRACT 

The simulated matter distribution on large scales is studied using core-sampling, clus- 
ter analysis, inertia tensor analysis, and minimal spanning tree techniques. Seven sim- 
ulations in large boxes for five cosmological models with COBE normalized CDM-like 
power spectra are studied. The wall-like Super Large Scale Structure with param- 
eters similar to the observed one is found for the OCDM and ACDM models with 
fi^/i = 0.3 & 0.245. In these simulations, the rich structure elements with a typical 
value for the largest extension of ^^(30 — 50)/i~^Mpc incorporate ~40% of matter with 
overdensity of about 10 above the mean. These rich elements are formed due to the 
anisotropic nonlinear compression of sheets with original size of ~(15 — 25)/i^^Mpc. 
They surround low-density regions with a typical diameter ~(50 — 70)h~^ls/lpc. 

The statistical characteristics of these structures are found to be approximately 
consistent with observations and theoretical expectations. The cosmological models 
with higher matter density O^ = 1 in CDM with Harrison-Zeldovich or tilted power 
spectra cannot reproduce the characteristics of the observed galaxy distribution due 
to the very strong disruption of the rich structure elements. Another model with 
a broken scale invariant initial power spectrum (BCDM) shows not enough matter 
concentration in the rich structure elements. 
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1 INTRODUCTION 

The phenomenon of Super Large Scale Structure (SLSS) was 
first observed as rare peculiarities in the visible galaxy distri- 
bution, with extreme parameters. Examples of this include 
the Great Void (Kirshner et al. 1983), the Great Attrac- 
tor (Dressier et al. 1987), the Great Wall (de Lapparent, 
Geller, & Huchra 1988; Ramella, Geller, & Huchra 1992), 
and the Pisces-Perseus supercluster (Giovanelli & Haynes 
1993). Several nearby superclusters of galaxies were de- 
scribed by Oort (1983a, b). SLSS was also found in deep pen- 
cil beam redshift surveys (Broadhurst et al. 1990; Willmer 
et al. 1994; Buryak et al. 1994; Bellanger & de Lapparent 
1995; Cohen et al. 1996) as rich galaxy clumps with typi- 
cal separations in the range of (60 — 120)/i~^Mpc {h is the 
Hubble parameter in units of 100 km s~^ Mpc~^). 

Recently the analyses of rich galaxy surveys with 
an effective depth ~(200 — 400)/i~-^Mpc, such as the 
Durham/UKST Galaxy Redshift Survey (Ratcliffe et al. 
1996) and the Las Campanas Redshift Survey (LCRS) 
(Shectman et al. 1996), have established the existence of 
the wall-like SLSS as typical phenomenon in the visible 
galaxy distribution incorporating ~(40 — 50)% of galaxies 



(Doroshkevich et al. 1996, hereafter LCRSl; Doroshkevich 
et al. 1997b, hereafter LCRS2; 1998a; 1999). The wall-like 
SLSS consists of structure elements with a typical diameter 
~ (30— 50)/i^^Mpc surrounding low-density regions (LDR) 

with a similar typical diameter Dldr ~ (50 70)/i~^Mpc. 

Within the wall-like SLSS elements the observed galaxy dis- 
tribution is also inhomogeneous (see, e.g., Fig. 5 of Ramella 
et al. 1992), and galaxies are concentrated in high density 
clumps and filaments. 

In LDR the galaxies are found to be concentrated within 
a random network of filaments. In distinction to the typi- 
cal wall-like superclusters, the galaxy distribution in LDR 
is predominantly 1-dimensional with a mean separation of 
filaments ~(10- 15)/i"^Mpc (LCRSl), and we caU this net- 
work Large Scale Structure (LSS). The LSS incorporates 
also ~50% of galaxies and is clearly seen in many redshift 
catalogues of galaxies (see, e.g., de Lapparent, Geller & 
Huchra 1988). These results extend the range of investigated 
scales in the galaxy distribution up to ~ 100/i~^Mpc. 

While the LSS was predicted by the nonlinear theory 
of gravitational instability (Zel'dovich 1970), and it was re- 
produced in simulations before its discovery in observations 
(see, for reference, Shandarin & Zel'dovich 1989), the ob- 
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servation of the rich and typically wall-like SLSS was quite 
unexpected. Thus, in simulations the representative SLSS 
was found only recently for a CDM model with low den- 
sity and a cosmological constant (ACDM, Cole et al. 1997). 
This simulation demonstrates that for suitable parameters 
of cosmological models the formation of the wall-like SLSS 
is compatible with the standard CDM power spectrum of 
Gaussian initial perturbations. 

The formation and evolution of structure on all 
scales are described by an approximate theoretical model 
(Demianski & Doroshkevich 1998a, b) based on Zel'dovich 
nonlinear theory of gravitational instability (Zel'dovich 
1970; 1978; Shandarin & Zel'dovich 1989). The model shows 
that the SLSS formation can be related to matter infall into 
the large wells of the gravitational potential of the initial 
perturbations. This model connects the structure parame- 
ters with the main parameters of the underlying cosmologi- 
cal scenario and the initial power spectrum. It shows that the 
impact of large scale perturbations is important throughout 
all evolutionary stages. In particular, the influence of these 
perturbations modulates the merging of smaller structure 
elements promoting the evolution within the SLSS elements 
and depressing it in LDR. 

The simulations are able to take into account this inter- 
action throughout all evolutionary stages and, therefore, are 
the most suitable way to study the LSS and SLSS properties 
and evolution. This implies, however, that the simulations 
need to be performed in very large boxes both to provide 
us with reasonable statistics of the walls and, of particular 
importance, to correctly describe the large scale part of the 
initial power spectrum of perturbations, responsible for the 
wall formation and the mutual interaction of small and large 
scale perturbations. In practice, this means that we need ex- 
treme parameters in the simulations. On the other hand, in 
order to compare various cosmological models, a broad set 
of simulations has to be prepared. 

Here results are presented from the analysis of simula- 
tions using five cosmological models. It is shown that the 
models with high matter density, O^ ~ 1, cannot repro- 
duce adequately the observed properties of the large scale 
matter distribution. The models with lower matter density, 
in particular the OCDM and ACDM models, seem to be 
more promising, since they can reproduce the general ob- 
served characteristics of both the LSS and the SLSS. In these 
models the rich structure elements, formed by a nonlinear 
matter compression, contain a significant matter fraction, 
frse ~ 0.4 — 0.5, which can be easily discriminated. A more 
detailed investigation of the nonlinear matter evolution on 
large scales is a further goal of this paper. 

The small scale matter clustering resulting in the de- 
struction of structure elements restricts the class of cosmo- 
logical models which are capable of reproducing the observed 
LSS and SLSS. The instability of a sheetlike matter dis- 
tribution similar to the observed and simulated SLSS was 
considered (in the linear approximation) by Doroshkevich 
(1980) and Vishniac (1983), and it was recently simulated 
by Valinia et al. (1997). 

Following our previous papers (LCRSl; LCRS2; 
Doroshkevich et al. 1997) we concentrate on the geometrical 
properties of the matter distribution, in particular, to the 
proper sizes and spatial distribution of filaments and wall- 
like structure elements. The popular correlation analysis is 



not so useful at scales > 10/i^^Mpc discussed below, and 
other techniques provide us with more essential results. We 
employ the core-sampling approach introduced by Buryak et 
al. (1994), the standard cluster analysis supplemented by the 
inertia tensor technique (Vishniac 1986; Babul & Starkman 
1992), the analysis of the variations of number of clusters vs. 
linking length (NCLL method, Doroshkevich et al. 1997b), 
and the minimal spanning tree (MST) technique (Barrow et 
al. 1985; van de Weygaert 1991). These methods were uti- 
lized recently for the investigation of structures in the LCRS 
(cp. LCRSl and LCRS2) and Durham/UKST redshift sur- 
veys (Doroshkevich et al. 1999). These results shall be used 
for the comparison with the structure parameters derived 
from simulations. The different methods are complementary 
to each other, and, thus, they allow us to characterize the 
simulated matter distribution in different important aspects. 

The observed distribution of galaxies and the simulated 
distribution of the DM cannot be identical as the galaxy for- 
mation is sensitive to additional factors (e.g., to the reheat- 
ing) and, moreover, galaxies, probably, mark only the high- 
est peaks of density perturbations. This means that some 
parameters of DM structure elements such as their overden- 
sity and proper sizes can differ from that found in observa- 
tional catalogues. The comparative analysis performed for 
one simulation (Doroshkevich et al. 1998a) confirms, that in 
some respect the spatial distributions of DM and 'galaxies' 
are different. A more detailed comparison of the observed 
and simulated matter distribution implies an identification 
of 'galaxies' in the simulated DM distribution. This means 
that a certain bias model needs to be utilized (see, e.g., dis- 
cussions in Sahni & Coles 1995, and Cole et al. 1998). Both 
problems are, however, equally important, and the distribu- 
tions of both the galaxies and the DM must be studied. 

This paper is organized as follows. The simulations and 
the analysis techniques utilized are shortly described in Sees. 
2 and 3. In Sec. 4 the general characteristics of the considered 
simulations are discussed that allow us to roughly discrimi- 
nate the cosmological models and to select the most realistic 
ones for a more detailed investigation. Our main results are 
presented in Sees. 5 and 6. Sec. 7 is devoted to the compar- 
ison with theoretical expectations of the DM distribution, 
and in Sec. 8 we discuss mock galaxy catalogues using some 
simple bias models. The conclusion and a discussion can be 
found in Sec. 9. 



2 SIMULATIONS 

We used five cosmological models as a basis for our anal- 
ysis. The COBE normalized SCDM model is taken as a 
reference model despite its difficulties in explaining already 
standard measures of galaxy clustering as the power spec- 
trum and the correlation function of galaxies and galaxy 
clusters (cp. e.g. Ostriker, 1993). Alternative models with 
Q.m = 1 include modifications of the primordial power spec- 
trum, in particular by introducing a tilt oc fc"'^ of the power 
spectrum (TCDM), or a break at a certain scale, (BCDM). 
Both are infiation motivated, using either an exponential 
inflation potential (Lucchin & Matarrese 1985), or a double 
inflation scenario (Gottlober, Miiller, & Starobinsky 1991). 
The BCDM is specified by two parameters, the location of 
the break at k^^^^y^ — 1.5/i^^Mpc, and the relative power 
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Figure 1. Simulated point distribution in rcdshift space for 
ACDM model at redshift z=0 



Table 1. Parameters of simulations. 



Nr, 



500 300=* 600^ 
400 2563 512^ 



h 

0.5 
0.5 

0.5 500 

0.5 500 3003 goO^ 0.60 714 
0.6 500 300^ 600^ 0.74 550 
0.7 500 300^ 600^ 1.12 913 

r3 c:i93 



300^ 



600^ 



0"8 O-Dei 

1.05 1157 

1.12 1127 

1.25 1293 



model Vim 

SCDMl 1 

SCDM2 1 

TCDM 1 

BCDM 1 

OCDM 0.5 

ACDMl 0.35 

ACDM2 0.35 0.7 400 256^ 512^ 1.30 873 
f2m is the matter density, h is the dimensionless Hubble constant, 
Lbox is the box size in /i^^Mpc, Np and N^^u are the particle and 
cell numbers, its is the mass variance for the linear spectrum at 
the scale 8/i~^Mpc, <T„ei is the variance of linear particle velocity 
in km/s. 

on both sides of the break, A = 3. These parameters were 
originally chosen to obtain optimal linear fits to the vari- 
ous large-scale structure observations (Gottlober, Miicket, 
& Starobinsky 1994), and later tested against A^'-body sim- 
ulations (Amendola et al. 1995; Kates et al. 1995; Ghigna et 
al. 1996; Retzlaff et al. 1998). Both TCDM and BCDM mod- 
els seemed to be promising since they have reduced power at 
galactic scales with respect to the COBE normalized SCDM 
model. 

Furthermore, two models are studied which are based 
on the wide range of observations pointing to a lower mat- 
ter density in the universe. First, we study an OCDM model 
with fim = 0.5 violating the inflationary paradigm of a spa- 
tially flat universe. Second, we take a model with a cosmo- 
logical constant, which has Q.m = 0.35 and a vacuum en- 
ergy leading again to a spatially flat universe. This ACDM 
model bears some advantage in alleviating the tight age con- 
straint of the universe. For all models the standard para- 
meterization of the CDM transfer function of Bardeen et al. 
(1986, hereafter BBKS) was used. In Table 1 the main pa- 
rameters of simulations are listed, including the matter den- 
sity Q.m, the dimensionless Hubble constant h, the box size 
Lhox, the particle number A'p and the cell number Nceii- Two 
models (SCDM and ACDM) were simulated with different 
resolutions and with slightly different initial amplitudes. 



The power spectra were normalized according to the 
two year COBE measurement following Bunn et al. (1995), 
(SCDMl, TCDM, BCDM, and ACDMl), or to the four 
year data according to the description of Bunn and White 
(1997), (SCDM2, OCDM, ACDM2), in both cases assum- 
ing pure adiabatic perturbations and a baryon content of 
Q.bh^ — 2 • 10^^ as predicted by big bang nucleosynthesis 
(see, e.g., Schramm 1998). Our later discussion will show 
that the smaller amplitude of the four year normalization, 
or a contribution of gravitational waves and/or other in- 
homogeneities to the COBE signal do not influence signifl- 
cantly the main conclusions. The amplitude of perturbations 
is characterized by the mass variance for the linear spectrum 
at the scale 8/i~^Mpc, ag,, and the 3-dimensional velocity 
dispersion, o^ei , gained from all dark matter particles in the 
simulation. 

The simulations were run in boxes of comoving size 
Lbox = 500/i~^Mpc and Li,ox = 400/i~^Mpc, respectively, 
to provide a good statistics of perturbations in the range 
of wave numbers k^^ ~ (10 — 30)ft^^Mpc responsible for 
the SLSS formation and, so, to improve the statistical char- 
acteristics of SLSS elements and the description of the in- 
teraction of small and large scale perturbations. For such 
boxes the formation of the majority of structure elements 
is described by higher harmonics of the primordial density 
waves, with I > 8 — 10. We use the PM code, described in 
more detail in Kates et al. (1995) and Retzlaff et al. (1998), 
with Np = 300^ or 256^ particles in Nceii = 600^ or5123 grid 
cells, respectively. These parameters provide a resolution 
~ 0.9/i"^Mpc and the mass resolution ~(1 - 3) • lO^M©. 

Most statistics can be calculated only for subsamples of 
the huge data sets. Therefore, we mostly used slices of thick- 
ness 50/i^^Mpc of the simulation box, that are about 10% 
of the complete volume. Even this volume provides us with 
a reasonable representation of the SLSS elements. This high 
stability of structure parameters is a direct consequence of 
the large box size used. To test the reliability, we repeated 
the analysis for different slices and constructed some statis- 
tics for the full volume taking the ACDMl model. The main 
difference between the analysis of the full sample and slices 
is a variations of the mean velocity dispersion of clusters by 
5—7 percent. Similar differences are characteristic of the dif- 
ferent realizations of the same cosmological model (SCDMl 
vs. SCDM2, ACDMl vs. ACDM2). Below we give the basic 
results for the larger simulations only, they are denoted as 
SCDM and ACDM. 

The comparison of simulated and observed parameters 
of the SLSS has to be done in the redshift space while the 
theoretical predictions are usually made in the comoving 
space. Hence, our analysis was performed twice. In the case 
of redshift space, we added an apparent displacement corre- 
sponding to the peculiar velocity of the particles along one 
axis divided by the Hubble constant. The comparison of re- 
sults from real and redshift space allows us to establish the 
influence of the velocity dispersion on the final estimates. 

Fig. 1 shows a wedge diagram of the simulated mat- 
ter distribution in redshift space for the ACDM model at 
redshifts z = 0. To each particle we assigned a luminos- 
ity chosen at random from a Schechter luminosity func- 
tion typical for the LCRS galaxies (Lin et al. 1996). We 
projected the particles into a wedge of angular extension 
80° X 1.5°, and we kept simulated galaxies in the magnitude 
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range 14 < m < 18, about the range of the LCRS. The 
mock sample contains about 7500 galaxies. The figure illus- 
trates that the overdense regions form almost a percolating 
system, with huge structured systems between radial veloc- 
ities of (10000 — 40000) km/s. A sparse filamentary matter 
distribution occupies the low-density regions. 



3 CORE-SAMPLING, CLUSTER ANALYSIS, 
AND MINIMAL SPANNING TREE 
TECHNIQUES 

In this paper we are interested in the investigation of large 
scale structure, and specific methods are to be used in or- 
der to characterize the simulated matter distribution. These 
methods are not so popular as, for example, the correlation 
analysis, but they are well defined and allow us to character- 
ize the matter distribution on large scales comparable with 
observations and with the scales predicted by theoretical 
considerations. 



3.1 Core-sampling approach 

The core-sampling method proposed by Buryak et al. (1994) 
was described in detail in LCRSl and tested with Voronoi 
tesselations by Doroshkevich, Gottlober, & Madsen (1997). 
It allows us to discriminate the filamentary and the sheet- 
like structure elements and to find two quantitative charac- 
teristics of the structure, namely, the surface density of fila- 
ments, (T/, that is the mean number of filaments crossing a 
randomly oriented unit area (i.e. 1 /i~^Mpc^), and the linear 
density of sheets, (Js, that is the number of sheets crossing 
the unit length (i.e. 1 fe~^Mpc) of a random straight line. 
These parameters are equivalent to the mean separation be- 
tween sheet-like structure elements, Ds, and filaments, Df. 

-1/2 



Ds 



l/(Ts 



D 



1 



'f 



(3.1) 



i.e. these lengths represent the mean free path between 
sheet-like and filamentary structure elements. 

The core-sampling method also allows us to determine 
the masses and velocities of structure elements that inter- 
sect a sampling core, i.e. it provides the surface mass density 
and dynamical characteristics of structure elements. These 
parameters are used for the comparison with theoretical ex- 
pectations. 

3.2 Minimal Spanning Tree technique 

The MST is a unique network associated with a given point 
sample and connects all points of the sample to a tree in a 
special and unique manner. Some definitions and capabilities 
of this approach are described by Barrow, Bhavsar, & Son- 
oda (1985) and van de Weygaert (1991). Here we will restrict 
our investigation to the analysis of the frequency distribution 
of the MST edge lengths Wmst{1) (the FDMST method). 
The potential of the MST approach is not exhausted by this 
application. It allows us to characterize, in particular, the 
morphology of structure elements and the typical size of the 
structure network. 

At large distances any correlations between the particle 
positions are small, and it can be expected, that the edge 
length distribution Wmst{1) is similar to that of a Poisson 



distribution. But for filaments this distribution will be dom- 
inated by a Poisson distribution with 1-dimensional support 
(ID), whereas for sheet-like elements a nearly 2-dimensionaI 
(2D) random point distribution is typical. This means that 
the function Wmst{1) can be used to characterize statisti- 
cally the dominant point distribution in the sample. To do 
this, the FDMST can be fitted to the six-parameter function 



WMsrix) 



-W. f^ e 
ax 



Fj(^) 



1/{Imst), (3.2) 



Ffix) ^ iPix"' + fSiX^'Y 



^^'^' ^ Ff dx ^ ^^ PixPi + I32XP2 ' 

where I and (Imst) are the edge lengths and the mean edge 
length of the tree, and Wo provides the normalization of the 
FDMST. The function F represents a power law both for 
small and large x, but it allows a continuous variation of 
the power index p{x) with the edge lengths x. 

Here we are mainly interested in the power index p{x) 
for larger x that characterizes the underlying geometry of 
the point distribution on large scales. A Poisson point dis- 
tribution with ID and 2D support is characterized by the 
power indices p — 1 and 2, respectively. Therefore, the 
asymptote of the function p{x) at large x characterizes the 
geometry of the structure elements, and it can be compared 
with similar parameters recently found for the observed 
galaxy distribution (LCRS2). This approach was tested with 
the simulations of ID, 2D and 3D Poissonian-liko point dis- 
tributions. 



3.3 Cluster analysis and variations of the number 
of clusters vs. linking length — the NCLL 
method 

The standard cluster analysis (friend-of-friends method) is 
used widely in numerical simulations and is well known (see, 
e.g., Sahni & Coles 1995). Here we employ it, first of all, in 
order to define the structure elements for a more detailed 
investigation of their properties. The cluster analysis can 
also be used for the description of the matter distribution 
on large scales. The function Wmst{1) is closely connected 
with the number of clusters N^i {ru„k) because 



tW( 



(3.3) 



Kl'irunk) ^Np WMST{l)dl, 

where A^p is the number of points in the sample under in- 
vestigation. Therefore, the FDMST and NCLL approaches 
are similar in many respects. 

We use a five-parameter fit of the cluster number vs. 
linking length relations, NJ:p (ru^k)- 



Nrj{b) 



Nr, e 



-Ff(b) 



/47r 3 \ 



1/3 



(3.4) 



Here Ff{b) is defined by (3.2), b is the dimensionless linking 
length, I3i,[32,pi,p2 and pg are dimensionless fit parameters, 
and Tip is the 3D number density of points. 

Here we are mainly interested in the variation of the 
power index p(6) at large b. This method is complementary 
to the FDMST analysis, and it allows us to get an indepen- 
dent fit to the power index p{b) at large b. 
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The NCLL method can also be extended, and, for the 
more detailed characteristics of the matter distribution, the 
variation of single particles, doublets, triplets and other clus- 
ters vs. linking length can be considered as well. In this pa- 
per, we consider only the total number of clusters N^l (rii„k) 
for the comparison with results obtained with the FDMST. 



3.4 Inertia tensor method 

The sizes of structure elements can be found with the inertia 
tensor method (Vishniac 1986; Babul & Starkman 1992). For 
each cluster the inertia tensor lij is expressed as 



N„ 



E 



JO) 



(a.-^Dfe-gr) 



(3.5) 



JO) 



where qt and q^ are the coordinates of the particles and of 
the center of mass of the cluster, respectively, Nmem is the 
number of cluster members. The conventional normalization 
5 has been taken to be consistent with a homogeneous el- 
lipsoid. Hence, the principal values of the tensor lij , namely 
the length (diameter), L, the width, w, and the thickness, 
t, {L > w > t), give us objective estimates of the cluster 
size and of the volume, V, and of the mean overdensity of 
cluster, S 



V- 



6 



-Lwt, 



Nrr 



(3.6) 



where Nmem is a number of point in cluster. 

This raw estimate is clearly of limited accuracy, but it 
is easy to calculate. The reliability of this estimate is high 
for compact regular clusters, and in general it depends on 
the cluster shape, in particular on its lumpiness (see, e.g., 
Sathyaprakash et al. 1998). We found from the simulations, 
that for large linking lengths, clusters are very lumpy in its 
outer regions. Then the ellipsoidal approximation leads to 
an artificial growth of the width and the thickness of the 
clusters. However, the cluster diameter L provides a stable 
characteristic of the cluster size. 



4 RICH STRUCTURE ELEMENTS IN 
SIMULATIONS 

The cluster analysis shows that rich structure elements 
(RSE) are usually represented by rather compact wall-like 
objects, and our methods give us more reliable information 
about their properties, some of which can be directly con- 
nected with the parameters of the cosmological model (see 
discussion in Sec. 7). In contrast, the discrimination and 
identification of poor structure elements is always difficult, 
as in the observed galaxy distribution they usually represent 
a filamentary component in a random network. This means 
that the discrimination and statistical description of such el- 
ements is often uncertain as their shapes are entangled due 
to many irregular branches. 

Because of this, in this paper we give the main attention 
to the RSE. Some statistical parameters of the filamentary 
component have nonetheless be found. They are discussed 
below. 



4.1 Identification of structure elements in 
simulations 

The large size of the computation box allows us to obtain a 
representative set of large clusters which can be associated 
with the observed RSE. The clusters were found for differ- 
ent richness thresholds Nthr, and for varying linking lengths 
riink- The linking length is directly connected with the den- 
sity threshold bounding structure elements, nthr- Indeed, 
no particles of a cluster are separated from the neighbor by 
more than the distance rank, therefore a lower limit to the 
cluster density is 



nthr > npb 



(4.1) 



where b is given by (3.4). 

For more interesting models the parameters Nthr and b, 
used for the identification of structure elements, as well as 
the number of identified structure elements, Nrse, are listed 
in Table 2. We also give the value bperc which corresponds to 
the linking length when the largest cluster of the sample ac- 
cumulates ~(25 — 30)% of the points. This is similar to, but 
not the exact percolation threshold (Klypin & Shandarin 
1992). The smooth shape of FDMST (no break or cutoffs, 
see the discussion in Sec. 6.6) shows, that in the simulated 
matter distribution even for large b, a significant fraction of 
points is not accumulated by the largest cluster. Neverthe- 
less, the fast growth of the largest cluster distorts the cluster 
properties, therefore smaller linking lengths must be used to 
obtain the typical characteristics of rich structure elements. 

A few parameters for the RSE are applied to discrimi- 
nate cosmological models, to select the most realistic mod- 
els for a detailed analysis, and to find a reasonable range of 
the threshold for the discrimination between RSE and LDR. 
These parameters, also listed in Table 2, are the fraction of 
matter accumulated by such structure elements, frse, the 
overdensity Srse given by (3.6), the velocity dispersion of 
these structure elements, ct„, and the velocity dispersion of 
matter within the structure elements along the three princi- 
pal axes of their inertia tensor, ai, a2, 0-3. The overdensity, 
the velocity dispersion of the clusters, and the inner veloc- 
ity dispersion were averaged over all clusters in the sample 
within the range of richness under consideration. Each clus- 
ter was weighted by the number of cluster members, Nmem, 
as this provides parameters typical for the majority of con- 
sidered points. These characteristics are sensitive to Qm and 
h, and to the amplitude and the shape of power spectrum. 
They are found for all simulations. The overdensity and the 
component of the inner velocity dispersion along the axis 
with the smallest component of the inertia tensor, 0-3, are 
plotted in Figs. 2 and 3 for SCDM, OCDM and ACDM mod- 
els. 



4.2 Overdensity and velocity dispersions of rich 
structure elements 

The simulated characteristics of the RSE must be com- 
pared with the observed parameters. The velocity dispersion 
within wall-like structure elements was roughly estimated by 
Oort (1983a) to be af" ~ (350 - 400) km/s. An estimate 
of the bulk velocity a^^" ~ 400 km/s can be taken from 
Dekel (1997). The fraction of galaxies accumulated by the 
RSE, frse, and the corresponding overdensity Srse, were es- 
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Table 2. Parameters of rich structure elements in comoving (cm) and rcdshift (rs) space. 



model 



^link 



b^ bL 



Nrs 



fr 



0-3 







/i-iMpc 












km/s 


km/s 


km/s 


km/s 




SCDM-cm 





0.75 


0.38 


2.22 


1370 


0.46 


43.1 


670 


690 


670 


590 




SCDM-rs 





0.75 


0.38 


1.95 


1134 


0.45 


10.2 


610 


716 


743 


703 




OCDM-cm 





0.95 


0.77 


1.39 


474 


0.39 


4.0 


416 


333 


341 


343 




OCDM-rs 





0.95 


0.77 


1.45 


469 


0.44 


3.8 


400 


325 


340 


344 




ACDM-cm 





1. 


0.90 


1.52 


752 


0.42 


38.1 


596 


610 


582 


514 




ACDM-rs 





1. 


0.90 


1.30 


697 


0.45 


14.1 


548 


565 


590 


530 




ACDM-cm 


1 


1. 


0.90 


1.37 


561 


0.24 


12.1 


496 


510 


527 


500 




ACDM-rs 


1 


1. 


0.90 


1.44 


599 


0.28 


10.1 


458 


473 


512 


497 




inkt b arc the threshold linking len 


5ths, bperc characterizes 


(approximate) percolation, Nrse 


frse and Srse 


arc the number, fraction of 



accumulated particles and mean overdensity, au is the dispersion of velocity of RSE and cri ,(T2 and 0-3 are the velocity dispersion of matter 
within the structure elements along the three principal axes of their inertia tensor. 
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Figure 2. Overdensity <5rse in redshift space vs. matter frac- 
tion concentrated within RSE, frse, for five richness thresholds: 
Nthr = 100 (dots), Nthr = 200 (dashed line), Nthr = 300 (stars), 
^thr = 500 (long-dashed line), Nthr = 1000 (open squares). Dot- 
ted lines show the observed parameters of RSE. 



Figure 3. Velocity dispersion 0-3 in redshift space vs. matter 
fraction concentrated within RSE, frse, for the same five richness 
thresholds. Dotted lines show the observed parameters of RSE. 



timated for the LCRS as frse ~ 0.4 — 0.5 and Srse ~ 10 
(LCRSl and LCRS2). The spatial distribution of RSE can 
be characterized by the mean separation of RSE along a ran- 
dom straight line, which was found as ~(50 — 60)^~ Mpc 
in the LCRS. Using these estimates, we can also find an 
approximate demarcation between RSE and LDR. 

The analysis shows that the mean velocity of the clus- 
ters (u) is negligible in comparison to the velocity dispersion 
au- This dispersion depends only weakly on the cluster rich- 
ness, but it is sensitive to parameters of the cosmological 
model. Because of the vortex-free character of the initial 
velocity field the velocity u measures the random variation 
of the gravitational potential over the cluster and, there- 
fore, the value (t„ is approximately proportional to the am- 
plitude of initial perturbations given by the mean velocity 
dispersion, a^ei, listed in Table 1. Theoretical considerations 
(Demiahski & Doroshkevich 1999b) describe this connection 
quantitatively (see also the discussion in Sec. 7) . It can be ex- 



pected that (Tu is slightly smaller in the MDM model where 
the fraction of hot DM particles makes the potential distri- 
bution more smooth. For the SCDM and TCDM models, the 
dispersions (t„ exceed the observed value by about a factor 
of 2, and by a factor of about 1.5 for ACDM model. For 
the OCDM and BCDM models, the simulated and observed 
velocity dispersions are in general consistent. 

The mean overdensity and the inner velocity dispersion, 
as, along the shorter cluster axis are very sensitive both to 
the cosmological model and to the subsample of clusters un- 
der investigation. For the most realistic models, they allow 
us to estimate a suitable range for the linking length, runk, 
and the threshold richness of clusters, Nthr- For three mod- 
els the overdensity, Srse, and the inner velocity dispersion, 
(73, are plotted in Figs. 2 and 3 vs. the matter fraction of 
the clusters, using five richness thresholds, Nthr- In order to 
compare the parameters of the clusters with observations, 
the analysis has been performed in redshift space. The ob- 
served estimates of overdensity, galaxy concentration and 
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velocity dispersion along the shorter principal axis of RSE, 
0-3, are plotted in Figs. 2 and 3 as well. 

The simulated values of the velocity dispersion (73 ex- 
ceed the theoretical expectations (Demiariski & Doroshke- 
vich, 1999b) due to the wall disruption and the formation of 
high density clumps. This is clearly seen from the isotropy of 
dispersion listed in Table 2 (we have always cri ~ (T2 ~ (T3). 
The rate of this disruption depends on the degree of matter 
compression and, therefore, on the mean overdensity of the 
clusters. 

The simulations with D,m = 1 (SCDM, TCDM and 
BCDM models) cannot reproduce the main observed char- 
acteristics of the RSE over the range of considered richness 
thresholds 100 < Nthr < 1000 and over a reasonable range 
of linking lengths. For the SCDM and TCDM models, a rea- 
sonable matter concentration is connected with a very large 
velocity dispersion. Therefore, even moderate changes in the 
power spectrum normalization of these models cannot im- 
prove the parameters of the rich clusters. For the BCDM 
model, a reasonable velocity dispersion is accompanied by 
a very small matter concentration within the RSE. For this 
model a large bias between the galaxies and the dark matter 
distribution (a large factor ~ 2 — 3 is required) could, in prin- 
ciple, reduce this disagreement. The same effect can also be 
reached by an increase of the amplitude of the power spec- 
trum by a factor of about two in comparison to the used 
normalization. Thus, already the first step of our analysis 
shows that probably the models with Om ~ 1 cannot be 
considered as realistic. 

Models with a smaller matter density fim ~ 0.5 
(OCDM) and Qm = 0.35 (ACDM) show better results. For 
the ACDM simulation reasonable structure parameters are 
found for 6'^ — 0.90 and Nthr = 200. Some excess in the 
velocity dispersion (t„ and in cri, (J2, and 173 (a factor of 
about 1.5) points to the overevolution of this model. Better 
results can be reached by the variation of the DM composi- 
tion and/or for smaller values of the cosmological parame- 
ters hilm- For the simulated OCDM model, all parameters 
of RSE at b^ = 0.77 and Nthr = 200 are consistent to the 
observed ones in the range of our precision. A deficit in the 
overdensity {Srse ~ 4) is within the range of a possible large 
scale bias, i.e. a higher concentration of luminous matter 
(galaxies) in rich structure elements in comparison to the 
concentration of DM (see Sec. 8). 

The considered low density models are most promis- 
ing for a detailed investigation. In these models the COBE 
normalization is also consistent with the observed charac- 
teristics of rich clusters of galaxy (cp., e.g. Cole et al. 1997; 
Bahcall & Fan 1998). The curves in Figs. 2 and 3 allow us 
to establish a rough boundary between the RSE and LDR, 
both in terms of the variables Nthr and rii„k (or 6), and 
in terms of physical variables frse and Srse ~ {rirse) / (rip). 
Furthermore, detailed investigation of RSE and LDR allow 
us to test and to correct this demarcation. 

A significant redshift dependence of the main parame- 
ters of RSE is found for ACDM models at < z < 3. Thus, 
already at z = 1 a smaller matter concentration is found, 
and at z = 3, the matter fraction in RSE, frse ~ 0.03 — 0.05, 
is negligible. This means that for these models the RSE are 
sensitive indicators of the initial amplitude of perturbations. 
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Figure 4. Mean separation of the RSE, Ds, and the surface den- 
sity of filamentary component, af, vs. the matter fraction con- 
centrated within the structures, /, in redshift space, and at 2 = 0, 
for the OCDM (dots) and ACDM (stars) models. 



5 MEAN SEPARATION OF FILAMENTARY 
AND SHEET-LIKE COMPONENTS IN THE 
SCDM, OCDM AND ACDM MODELS 

In this Section, properties of simulated structures are ex- 
amined with the core-sampling method. This method allows 
us to find the mean separations between filamentary and 
sheet-like structure elements, respectively. The analysis of 
samples obtained by systematic rejection of sparser struc- 
ture elements allows us to reveal the characteristics of typ- 
ical structures. These data can be compared with similar 
results obtained for the LCRS (cp. LCRSl). The analysis 
was performed for the SCDM and OCDM models at 2 = 0, 
and for the ACDM model at z = and z = 1. 

For the core-sampling analysis a set of 196 cylinders 
with a radius of 1.7fe~^Mpc was prepared both in comoving 
and in redshift space. The mean number of points within the 
cylindrical cores amounts ~(400 — 600). The analysis was 
performed for 16 values of the cylinder radius, 1.7/i^^Mpc 
> Tcyi > 0.7/i^^Mpc. The separation of sheet-like elements, 
Ds, and the surface density of the filamentary component 
af — DJ^ , are plotted in Fig. 4 versus the fraction of matter 
/, remaining after rejection of sparse structure elements. 
It shows the OCDM and ACDM models in redshift space 
at z = 0. The parameters typical for the reliable structure 
elements are marked in Fig. 4 by dotted lines and listed 
in Table 3 together with similar parameters obtained for 
the LCRS (LCRS-1). Results are found to be close both 
in comoving and redshift space, and they coincide with the 
parameter range estimated for the LCRS. 

In all cases there is a clear signal from the SLSS com- 
ponent, but in contrast with results found for the LCRS, 
Ds increases slowly for small / < (0.6 — 0.7). This effect 
is probably caused by the variation of the covering factor 
of the sheet-like component, as discussed by Ramella et al. 
(1992) and Buryak et al. (1994). This effect is less promi- 
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Table 3. Structure parameters with the core-samphng method 
in redshift space. 



model z 

LCRS 0. 

SCDM 0. 

OCDM 0. 

ACDM 0. 

ACDM 1. 



JVy 



Ds 



Df 






thr l^s Uf 

/i-^Mpc h-^Mpc /i^iMpc 

55. ± 7 26.7 ±3. 10.5 ± 1. 

8 38. ± 9 19.2 ±3. 7.5 ± 2. 

6 48. ±11 19.2 ±3. 10.0 ±3. 

4 63. ±15 26.7 ±5. 12.0 ±4. 

4 45. ±11 23.6 ±3. 7.8 ± 2. 



The threshold A^t^^ is the minimal richness of the considered 
structure elements, Ds, Df, and D'T'''" are the mean separation 
of sheet-like elements and rich and poor filaments. 



nent for the OCDM model for which the disruption of RSE 
measured by the velocity dispersions, cri, (72, and (73, is also 
less significant. The weak variation of Ds{f) and the quick 
drop of af with / shows that a significant matter fraction 
(~ 0.4 — 0.5) is associated with the high dense sheet-like 
component that is also consistent with the observational re- 
sults in the LCRS (LCRSl). A population of rich filaments 
can be identified in Fig. 4, but in contrast with results from 
the LCRS, it is less representative. For the SCDM model 
we also found a strong variation of the number of particles 
from core to core. This is probably caused by the strong 
disruption of structure elements. 

The parameters D™'" listed in Table 3 correspond to 
the minimal separation of filaments found in our analysis at 
/ ^ L As it is common in simulations, the identification of 
poor filaments is difficult, and their parameters depend on 
the rejected background. It is well known that for a CDM 
power spectrum, very small DM pancakes form early, and 
the estimates of the minimum pancake size in simulations 
reflect the resolution (in our case given by the size of the 
computation cells). This means that D™™ listed in Table 3 
generally characterizes the procedure of background rejec- 
tion rather than properties of poor structure elements. This 
problem was discussed in detail by Doroshkevich et al. (1997 
& 1998b). The uncertainties in our estimates of Ds and Df 
are larger than those obtained for the LCRS. This shows 
that simulated DM structures are not fully consistent with 
the geometrical model on which the core-sampling method 
is based. 

The comparison of structure parameters found for the 
ACDM models at 2: = and z — 1 shows the moderate evolu- 
tion of structure which is approximately consistent with the- 
oretical expectations (Demiaiiski & Doroshkevich 1999b). 



6 PROPERTIES OF RICH STRUCTURE 
ELEMENTS 

The analysis performed in Sec. 4 and 5 shows that, in fact, 
for a significant matter fraction frse -^ 0.4, a strong nonlin- 
ear matter compression results in the formation of massive 
high density RSE. The existence of such RSE is a very es- 
sential feature of the large scale matter distribution. A more 
detailed analysis and the statistical description of RSE is de- 
scribed in this Section. 



6.1 Discrimination of rich structure elements 

The subsamplos of RSE were discriminated and analyzed at 
redshift z = in comoving and redshift space for the ACDM, 
OCDM and SCDM models, and for estimating the evolution, 
a similar subsample of RSE was also analyzed at redshift 
z = 1 for the ACDM model. As described in Sec. 4, the RSE 
were identified with rich clusters found for a suitable linking 
length, riink, and a richness threshold, Nmem > Nthr- The 
probability distribution functions discussed below depend 
on the used definition of structure element, and in our ap- 
proach, on these two parameters. The employed parameters 
of the subsamples are listed in Table 4. 

In the SCDM model, RSE are defined with linking 
lengths b'' --^ 0.48, that corresponds to a threshold density 
of clusters nthr > 2.1(np), whereas for OCDM and ACDM 
models b^ ~ 0.9 provides better results. For these samples, 
a significant fraction frse ~ 0.4 of all matter is contained 
in massive overdensity clumps with Nmem > Nthr = 200. 
The chosen values of 6^ are about 2 — 4 times smaller than 
the percolation threshold, bperci also listed in Table 2. This 
means that the RSE defined by these thresholds are actually 
isolated. Some parameters discussed below, in particular the 
comoving sizes, depend on the chosen b^ , and they increase 
for larger fe"^ and /or larger Nthr- 

These parameters are the basis for our selection of RSE 
in simulations. A more detailed comparison of observed and 
simulated RSE implies an identification of 'galaxies' in the 
simulated DM distribution (see Sec. 8). 

In the LCRS a similar fraction of galaxies, fgai ~ 0.4, 
is concentrated in clusters defined with a threshold linking 
length of b^'^ ~ (1 — 0.5), and a mean overdensity in the 
RSE of ~10 (LCRS2). This difference between the threshold 
densities of RSE used in simulations and in the LCRS is 
caused in main by the construction of the LCRS as a set of 
six slices with angular size 1°.5 and the effective thickness 
~(6 - 7)/i"^Mpc only. 

The impact of the slice thickness was tested using the 
mock catalogues prepared by Cole et al. (1998). It was found 
that the random intersections of RSE with relatively thin 
slices result in an artificial destruction of selected RSE. It is 
seen as the growth of the threshold linking length required 
for the selection of RSE and as stronger random variations of 
RSE properties with the linking length. Thus, it was found 
that for the full mock catalogues about (40 — 45)% of 'galax- 
ies' are incorporated in RSE already at 6^ ~ 1 that is close 
to the values used above. In contrast, for the slices with an- 
gular size 1°.5 the same 'galaxy' concentration in RSE is 
reached for b'^ ~ (1.5 — 2) that is comparable with that used 
for the LCRS (LCRS2). The small slice thickness depresses 
also the percolation and restricts the sizes of the richest RSE 
in the LCRS. The random mixture of fields observed with 
50 and 112 fibers (in four slices of the LCRS) amplifies this 
artificial destruction of selected RSE. 



6.2 Mass functions of rich structure elements 

The mass function of rich structure elements is written as 



Wm = 



Nr, 



-Nrse{l^), V = NraemKNraem) (6.1) 



where Nmem and A'j, are the numbers of points in the clus- 
ter and in the sample as a whole, Nrse is the number of 
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Figure 5. Mass function, IVm, for the OCDM, ACDM and 
SCDM models in redsliift space. Samples and models parameters 
are listed in Tables 1, 2 and 4. 
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Figure 6. The proper sizes of the RSE, X/(\),ui/(ui), and 0/{e), 
and the size of 'proto walls', Dprw/{Dprw), vs. the mass (richness) 
of the element, u = Nmem/ {Nmem) for the ACDM model in 
redshift space at z = 0. 



RSE elements for a given richness v, and (Nmem) is the 
mean number of points per RSE element. This function de- 
pends on the degree of matter concentration in RSE, and 
on the disruption of structures. This disruption divides the 
large structure elements into a system of high density clumps 
bridged by low density regions, and it also increases the frac- 
tion of low mass elements. The rate of disruption depends 
on the density contrast and, therefore, it is sensitive to the 
cosmological parameters. 

The mass functions are plotted in Fig. 5 for the OCDM, 
ACDM and SCDM models as determined in redshift space. 
The general character of the functions is very similar in these 
three models. Some excess of RSE in the tail, i.e. for masses 
f > (3-4), are found in the OCDM and ACDM models. The 
threshold richness cuts off the distribution at small masses, 
and it influences the value of the mean richness of the RSE, 
(Nmem)- The mean values are listed in Table 4 for all models 
under investigation. 

These mass functions are in general similar to the ob- 
served one in the LCRS (LCRS2), and at z^ < 6 they can 
also be fitted by a simple exponential law. In both cases, 
there are a few huge clusters with mass v ^ {5 — 10) which 
accumulate ~10% of points in the simulations. The rejection 
of low mass structure elements has a strong influence on the 
extent of the mass function. 



6.3 The proper sizes of the rich structure 
elements 

The proper sizes of rich structure elements are found with 
the inertia tensor method, applied to the subsample of RSE. 
All proper sizes depend on the mass of the RSE. The scaling 
can be approximated by 

,1/3/, 



L = u'/'X, 



t = v^ 



(6.2) 



Q.2<V = Nmem/ {Nmem) < 7 



The mass-averaged length (A), width (uj), and thickness (6) 
of clusters are listed in Table 4 together with the corre- 
sponding dispersions. Scatter plots of these values vs. the 
richness of RSE f for the ACDM model are shown in Fig. 
6. The distributions of proper sizes X, u;, and 9, are similar 
to Gaussian distributions with mean values and dispersions 
listed in Table 4. 

The scaling (6.2) describes well the mass dependence 
of the proper sizes which results from the relatively regu- 
lar shape of RSE and the moderate influence of boundary 
effects. In the LCRS a similar scaling is found for the two 
larger sizes, whereas the small size is weakly dependent on 
the richness (LCRS2) . This is probably caused by the special 
construction of the LCRS as a set of thin slices. The strong 
richness dependence found for the Durham/UKST redshift 
survey (Doroshkevich et al. 1999) could be caused by the 
relatively small size of the survey. 

Results listed in Table 4 show that for the OCDM and 
ACDM models, all proper sizes are found to be close (in 
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Table 4. Proper sizes of rich structure elements. 



model 

LCRS 
SCDM-cm 

SCDM-rs 
OCDM-cm 

OCDM-rs 
ACDM-cm 

ACDM-rs 
ACDM-cm 

ACDM-rs 



rink {Nmem) (A) a^ (uj) 



0.75 
0.75 
0.95 
0.95 
1.0 
1.0 
1.0 
1.0 



731 
870 
907 
998 
801 
919 
588 
641 



25.3 
9.2 
22.8 
24.2 
24.4 
14.2 
21.8 
19.6 
19.7 



6.0 
4.0 
6.6 
7.8 
6.8 
6.2 
5.5 
6.1 
4.4 



12.3 
4.3 
7.0 
12.1 
12.5 
6.2 
9.9 
7.8 
9.1 



3.1 

1.4 
2.6 
4.4 
3.2 
2.1 
2.9 
1.8 
2.1 



(e) 

5.7 
2.7 
3.6 
6.0 
7.2 
3.6 
5.0 
4.6 
4.8 



1.3 
0.6 
0.9 
1.9 
1.7 
0.9 
1.4 
1.3 
1.1 



\LJpru)/ 

26.4 
16.4 
16.0 
17.3 
18.3 
22.2 
23.2 
21.4 
21.1 



All mean sizes and dispersions are given by (6.2), (6.4) and Sec. 6.5 and are measured in h -"^Mpc 
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the range of the dispersion) to the sizes observed in the 
LCRS. Moderate variations of the threshold hnking lengths 
do not change the mean characteristics of the RSE, but the 
sizes of the largest structure elements are sensitive to such 
variations. 

For the SCDM model, all mean proper sizes in redshift 
space exceed the sizes found in comoving space. These differ- 
ences become smaller for the ACDM and the OCDM models, 
and they decrease with increasing redshift as it is shown by 
the comparison of the ACDM model at z = and z = 1. It 
is explained, in part, by the well known influence of the ve- 
locity dispersion ('finger of God effect'). But this effect also 
depresses the small scale clustering and, apparently, partly 
cancels the disruption of large structure elements typical for 
the nonlinear evolution of pancakes. This is clearly seen in 
the growth of the mean length (A), which exceeds any ve- 
locity dispersions and, therefore, is insensitive to the finger 
of God effect. In contrast, the influence of the small scale 
clustering on the thickness of structure elements is not so 
strong, and its growth in redshift space is probably caused 
by the direct influence of small scale velocities. Hence, we 
can take the value (L) in redshift space as a genuine diam- 
eter of the RSE, and the value (i) from comoving space as 
real thickness. For the middle size, (w), both effects can be 
important. 



6.4 The 'proto-size' of the rich structure elements 

The measured parameters of RSE allow us also to estimate 
the volume which initially contained all matter in RSE, and 
the degree of compression connected with its formation. A 
simple model was used for this estimate in LCRS2. As the 
lengths of the RSE exceed the two other sizes, we can neglect 
the matter compression along this axis, and we can consider 
the RSE formation as a 2D matter infall to the gravitational 
well. In this case we can obtain a simple estimate using the 
mass conservation: 

(nrsetw) Ri {np){toWo) « {np)Llr^ (6.3) 

J- 2 '^^^mem j- ^ j 

nL{np) 

where Urse and (rip) are the comoving number density of 
particles within the RSE and in the sample as a whole, re- 
spectively, t and w are the thickness and width of the RSE, 
and to, Wo and Lprw = V'^oto are the typical sizes of the 
'proto-structures', defined in the initially almost homoge- 
neous matter distribution. The mass dependence can be de- 



n-^. : 




Figure 7. The frequency distribution {Dsep)dW/dDsep of sep- 
aration of the RSE in comoving (panels a and c) and redshift 
(panels b and d) space at z=0 (panels a and b) and z=l (panels 
e and d). The exponential fits are shown by dotted lines. 



scribed similarly to Eq. (6.2): 



r _ 1/3 n 



(6.4) 



distributions with the mean values, (Dprw), and dispersions, 
aprw, listed in Table 4. The comparison of to and with 
Dprw shows that at z = the RSE formation can be roughly 
described as an asymmetric matter compression by a factor 
of ~2 along the middle axis («;), and a factor of ~(4 — 
5) along the shorter axis (t). At z — 1 the corresponding 
compression factors are ~1.5 times smaller. For the high 
density clumps at redshift 2 = 3 in the ACDM model the 
size of the 'proto-structures' is 



Lpr-u^iz = 3) ^ (18 ± 3.5)/i"^Mpc. 



(6.5) 



The values Dprw are plotted in Fig. 6 for given samples 
of RSE of the ACDM model. The parameters are very stable, 
in particular, they depend only weakly on the sample under 
investigation, and on the linking length, since usually the 
growth of the cluster sizes is accompanied by a drop in the 
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number density. But they are sensitive to the possible bias 
between the spatial distribution of DM and galaxies. 



6.5 The mean separation of the structure 
elements 

The mean separation of RSE can be found, applying a 
simple version of the core-sampling method for subsam- 
ples of RSE. A sample of 250 rectangular cores with sides 
10/i~^MpcxlO/i^^Mpc, oriented along one axis (the axis 
where the redshift distortions are added) and containing all 
particles of RSE, was prepared for the ACDM and OCDM 
models both in comoving and in redshift space for z — and 
for the ACDM model at 2 = 1. All particles are projected on 
the axis of the cores, and they are collected in clusters with 
the linking length, ri„k, used for the RSE preparation (Table 
4). The mean separation between these clusters (the 'mean 
free-path' between the RSE), (Dsep), and the dispersions, 
a Sep, are listed in Table 4. The large dispersion - more than 
50% of mean value - is typical for an exponential Poisson- 
like distribution. The frequency distribution of cluster sepa- 
rations is plotted in Fig. 7 together with the best exponential 
fit. The mean separation of the RSE depends weakly on the 
redshift, and it is consistent with the mean separation of 
wall-like elements found with the core-sampling approach 
in Sec. 5. For the OCDM model the numerical estimates 
are consistent to that found in the LCRS (cp. LCRS2), but 
for the ACDM model they exceed the observed values by a 
factor of ~ 1.5. 

The mean separation of RSE (or the ID number den- 
sity) allows us to obtain an independent estimate of the 
relative richness of RSE for a given subsample, i.e. for some 
linking length b and threshold Nthr, or a certain matter frac- 
tion in RSE, frse- The difference found above for the ACDM 
model indicates that in this case the same fraction of parti- 
cles is concentrated in a smaller number of RSE. To elimi- 
nate this difference, a threshold richness Nthr ~ 150 could 
be used, which has small effects on the other parameters of 
RSE. 

The mean separations between structure elements of the 
filamentary component and their dispersions, (Df) and ao, 
can be found through a similar analysis of subsamples pre- 
pared by removing all RSE from the full sample. We obtain 
(Df) ~ 9/i"^Mpc, CTD ~ 6.5/i"^Mpc, for the OCDM model, 
and (Df) ~ Uh'^Mpc, an ~ 10.5/i"^Mpc, for the ACDM 
model. This data agree well with the value D™'" listed in 
Table 3 for the LCRS. 



Table 5. Fit parameters for the FDMST and for the cluster distri- 
bution, pt, for full samples, RSE and LDR for the ACDM model. 



sample (Imst) ^mst Pmst 



Vt 





z = 


= 0, comoving space 




a TOT 


0.72 


0.72 0.60 ±0.02 


0.59 ±0.03 


bRSE 


0.30 


0.62 1.11 ±0.02 


1.02 ±0.05 


cLDR 


1.00 


0.80 0.83 ±0.02 


0.99 ±0.03 




z 


= 0, redshift space 




dTOT 


0.78 


0.78 0.91 ±0.03 


0.85 ± 0.04 


eRSE 


0.42 


0.83 1.60 ±0.04 


1.83 ±0.03 


f LDR 


1.00 


0.80 1.11 ±0.03 


1.11 ±0.02 




z = 


= 1, comoving space 




TOT 


0.90 


0.90 0.76 ±0.03 


0.97 ±0.03 


RSE 


0.35 


0.69 1.0 ±0.03 


1.30 ±0.05 


LDR 


1.00 


0.80 0.87 ±0.03 


1.10 ±0.02 




z 


= 1, redshift space 




TOT 


0.91 


0.91 0.97 ±0.02 


1.02 ±0.02 


RSE 


0.45 


0.89 1.4 ±0.05 


1.93 ±0.02 


LDR 


1.00 


0.80 1.1 ±0.03 


1.20 ±0.03 




1/<1„I 



Figure 8. The frequency distribution of edges of MST for the 
ACDM model. Sample are marked and sample parameters are 
listed in Table 5. The fit of FDMST with relation (3.1) is shown 
by solid line. 



6.6 Inner structure of RSE and LDR 

The point distribution in a sample can be characterized with 
the NCLL and the MST techniques as described in Sec. 3. 
These methods allow us to discriminate between the dom- 
inance of filamentary and sheet-like structures in a point 
sample, and to characterize the point distribution within 
separate structure elements. In this sense, the methods are 
complementary to the core-sampling approach. Thus, for the 
filamentary component, power indices pt and pmst ~ 1 can 
be expected, whereas for the sheet-like component, the ap- 
pearance of power indices pt and pmst ~ 2 seems to be more 
typical. 



Here we apply these methods to the ACDM model. We 
analyze separately the full sample, the RSE, and the LDR, 
the latter are obtained by removal of RSE from the full 
sample. The main results are presented in Figs. 8 and 9, 
and they are collected in Table 5 where {Imst) is the mean 
edge length of the MST and 



,47r 



hMST = {^^{np}) (Imst) 



(6.6) 



Here (up) is the mean number density in the sample under 
consideration. 

The FDMST are plotted in Fig. 8 for the full sample 
and for the RSE and LDR, both in comoving and in redshift 
spaces, and ai z = 0. The variation of the power indices vs. 
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the edge lengths are plotted in Fig. 9. The left-hand side of 
the FDMST describes the matter condensation within high 
density clumps that form the inner structure of filaments 
and walls. It is similar for all subsamples. The right-hand 
side of the FDMST characterizes the relative positions of 
these clumps and other particles of the subsamples. 

These figures demonstrate that in the comoving space 
for all samples, pmst ~ 1 for edge lengths 1/{Imst) > 0.5 — 
1. The NCLL method confirms that pt ~ 1 is reached for 
6 > 0.5, this is about 1.6 times smaller than b ~ 0.8 used 
for the RSE discrimination, and ~2.5 times smaller than 
the value feperc listed in Table 2. These values emphasize the 
joint character of the point distribution both within RSE 
and LDR which can be interpreted as a predominantly ID 
Poisson distribution typical for filaments. The mean edge 
length of the MST, {Imst), in the RSE is ~ 2 - 3 times 
smaller than within LDR but variations of bMST do not 
exceed ~20%. In redshift space the impact of the velocity 
dispersion erases the small scale structure of RSE. Therefore 
we find an apparent particle distribution similar to a 2D 
Poisson distribution. 

The characteristics obtained both with the NCLL and 
MST methods are similar to each other within the range of 
statistical uncertainty. In any case, independent on the geo- 
metrical interpretation, the power indices and typical scales, 
that characterize the spatial matter distribution in the RSE 
and LDR, are essentially different, and this verifies the ac- 
cepted discrimination of these regions. These results agree 
well with estimates of the mean overdensity, Srse, listed in 
Table 2 and obtained in other way. They confirm the essen- 
tial concentration of high density clumps in RSE. 

The results obtained for the LDR are consistent with 
the dominance of a filamentary component. For the RSE, 
the unexpected value of the power index in the comoving 
space can be considered as an indirect evidence in favor of 
RSE formation from earlier formed filaments. It can also be 
traced back to different factors, such as the wall disruption, 
the limited resolution of simulations and, therefore, further 
investigations are required. 

In the LCRS a power index pt ~ Pmst ~ 1.7 has been 
found for RSE, and pt ~ Pmst ~ 1 for LDR and for the total 
sample (LCRS2). This is comparable with our results for 
the redshift space listed in Table 5. The complicated inner 
structure of RSE is also seen in the LCRS and, more clearly, 
in the galaxy distribution within the Great Wall (Fig. 5 in 
Ramella, Geller, & Huchra 1992). 



7 CHARACTERISTICS OF THE EXPECTED 
AND SIMULATED DM STRUCTURE 

The reproduction of the main observed characteristics of 
the RSE in simulations with a standard CDM-like power 
spectrum verifies that the observed structure was formed 
during the nonlinear evolution of small initial perturba- 
tions, and, so, the characteristics of structure can be ex- 
pressed through the parameters of a suitable initial power 
spectrum of Gaussian fluctuations for a specific cosmolog- 
ical model. Statistical characteristics of the DM structure 
based on Zel'dovich nonlinear theory of gravitational insta- 
bility (Zel'dovich 1970, 1978; Shandarin & Zerdovichl989) 
were discussed by Domianski & Doroshkevich (1999a, b). 



1, h 

0, L_ 

2. h 

1, h 

0. U 

2, ^ 



0,0 0.5 1.0 1.5 2.0 

1/<1mst> 

Figure 9. Variations of power indices pmst for the same samples 
of ACDM model as in Fig. 8. Sample parameters are listed in 
Table 5. The fit of FDMST with relation (3.1) is shown by solid 
line. 



The approximate expressions, derived there, connect some 
of such characteristics to the parameters of the cosmologi- 
cal model. The comparison of the approximate analytic re- 
sults and the simulations reveals both the influence of factors 
omitted in the theoretical description and of random factors 
distorting the simulated structure. 

The surface density of RSE and velocity dispersions 
within RSE, and the velocity of structure elements seem 
to be most interesting. These values can be found using the 
simple version of core-sampling described in Sec. 6.5. We 
characterize the mass of each cluster within the rectangular 
sampling cores by the surface density of structure elements, 
m^, and the velocity dispersion within clusters also along 
the core, ar- These values can be found for the RSE. We con- 
sider also the dispersion of ID velocity of clusters along the 
core, aq, which can be found for the full samples, and for the 
RSE and the filamentary subsamples of structure elements 
separately. Some of these characteristics can be compared 
with similar 3D characteristics discussed above for RSE that 
demonstrates the influence of the sample selection and the 
averaging procedure. 

The theoretical parameters are expressed through the 
typical length scale, Iq, linked to the initial power spectrum, 
and typical dimensionless 'time', to, linked to the velocity 
dispersion, avei, and the parameter as listed in Table 1 as: 



^n 



To = 



kT{k)dk, lo 



6.6 , _i- . 
-—h ^Mpc, 



CTvel 



VspioHo ' 



/3«(0.43 + 0.57n„)-\ (7.1) 



TO = 0.55a8(f7™;i)"-^3« [l + 5.657(f7™/0'-']"'''' , 
where k is the comoving wave number, kmin 
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Figure 10. The distribution of mass surface density of RSE. Tlic 
best fit (7.3) is sliown by solid lines, the best fit by power laws is 
plotted by dotted lines. 



Table 6. Characteristics of mass and velocity distributions in 
comoving space for the full samples, RSE and LDR for the SCDM, 
OCDM and ACDM model. 



sample 



fir 



Tr 



SCDM 
TOT 501 230 

RSE 442 492 1.4 0.50 1.29 0.44 
LDR 508 151 - - 



0.66 
0.65 
0.66 



OCDM 

TOT 263 98 - - - - 0.30 

RSE 245 250 0.8 0.35 0.41 0.25 0.28 

LDR 277 76 - - - - 0.32 



ACDM 



TOT 360 155 
RSE 351 412 
LDR 374 119 



- - - - 0.42 
1.3 0.44 0.63 0.31 0.40 

- - - - 0.44 



the velocity dispersions (Tr and (Tq are given in km/s. 

2n/Lbox, kmax = kminN^^i^i, and T{k) is a transfer func- 
tion. The 'time' ro characterizes the ampUtude of perturba- 
tions and the reached period of structure evolution. For the 
SCDM, OCDM and ACDM simulations we have 



lo = 13.2, 22, 



26.9 /i^^Mpc, 



TO = T^ei = 0.43, 0.27 & 0.37 



TO = rg = 0.81, 0.31 & 0.46. 



(7.2) 



and differences between values r^ei & ts characterize the 
sensitivity of various estimates of amplitude to the small 
scale matter clustering. 

For Gaussian initial perturbations the distribution of 
pancake velocities along a core is expected also to be Gaus- 




2.5 



Figure 11. The distribution of velocity dispersions ar for the 
RSE. The best fit (7.3) is shown by solid lines. 



sian with a negligible mean velocity and a dispersion 

O-g ^ (T„ei/3 = HqIoPto/Vs. (7.3) 

Here a random orientation of pancakes with respect to the 
sampling cores is taken into account. 

In this case the distribution of surface density of RSE, 
m^,, can be expressed as follows: 

iV„ = ^e-^er/(Vi), x=^-P^, (7.4) 



fJ'VJ — 



bw{np)lo 



{rrin,) 



8(0.5 + l/7r)ro' = 6.6ro^ 



where a^ and b^ are fit parameters. The expected distribu- 
tion of the filamentary component is described by a more 
cumbersome relation. 

An expression as (7.4) with parameters a^ and br and 
X = brOr/ {(Jr) dcscribcs also the distribution of the ID inner 
velocity dispersion along the core, Or- For RSE this disper- 
sion can be expressed through Us and tq as 



Mr = 



{rrir) 3.3ro 



brHokli 2^ {np)/o x/3 



To 



3.3crq 

Holo/3 



(7.5) 



Using this relation we can measure the mean surface mass 
density of RSE (rrir) by the velocity dispersion. 

These relations allow us to estimate the model param- 
eters lo and To using the measured surface density of RSE 
{rriw}, the velocity dispersion within RSE, ar, and the veloc- 
ity dispersions of various populations of structure elements, 
Ug. To suppress the impact of small scale clustering, the 
analysis was performed for the subsample of RSE with the 
core side 10/i^^Mpc. Because of the small separation of fila- 
ments, for the full sample and subsample of LDR, discussed 
in Sec. 6.6, a core side 4/i~^Mpc was used. 

For RSE the measured distributions of surface density , 
rrin,, and velocity dispersion, ar, are plotted in Figs. 10 and 
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11 together with the best fits (7.4). The measured param- 
eters are Usted in Table 6, the theoretical expectations are 
given by Eq. (7.2). The distribution of pancake velocities is 
well fitted to the Gaussian function with dispersion aq listed 
in Table 6. 

These velocity dispersions are smaller by about 20% 
than those listed in Table 2 that demonstrates the influence 
of averaging procedures in getting the final parameters. The 
values To averaged over the measurements listed in Table 6, 
are 

(to) ~ 0.58 ± 0.09, ^ 0.29 ± 0.03, ^ 0.40 ± 0.05, (7.6) 

for the SCDM, OCDM, and ACDM models. Differences be- 
tween these values and estimates of tq given in (7.2) charac- 
terizes the influence of random factors such as the selection 
and disruption of structure elements and the real precision 
reached. 

The measured distributions of surface density of RSE 
are also well fltted by power laws with exponents k ~ 1.7 
that may be caused by the strong wall disruption. The power 
law can be reproduced analytically, assuming a set of clus- 
ters with spherically symmetric surface densities falling off 
according to a power law ads oc r^^. In this case the mass 
function of clusters in cores is also expressed by a power law 
as 



'els 



(7.7) 



with an exponent n — 2/7. For n ^ 1.7 we have 7 ^ 1.2 
that is close to the well known King's law, a^is oc (1 + 

,7 
elliptical galaxies 



r^ /r'^) ^ , 7 ~ 1, widely used to flt to the density profile of 



8 ANALYSIS OF MOCK CATALOGUES 

For a DM dominated universe the analysis of DM structures 
is very important in itself as the feedback of baryonic matter 
and galaxy formation to the DM evolution on scales larger 
than the mean intergalactic separation is small. On the other 
hand, almost all observed characteristics of the large scale 
structure are obtained for the galaxy distribution alone, and 
the galaxy distribution may be biased in comparison with 
the distribution of DM. Further, all observed galaxy cata- 
logues suffer from selection effects that influence our cosmo- 
logical conclusions. 

The selection effects are well studied and reproduced 
in available mock catalogues (see, e.g., Cole et al. 1998). A 
preliminary analysis of these catalogues reveals, for exam- 
ple, some distortions of observed parameters of RSE caused 
by the small angular size of the LCRS (see discussion in Sec. 
6). More detailed analysis of these catalogues with the tech- 
nique described above is in preparation. It allows us to find 
the optimal strategy of data analysis which suppresses the 
influence of selection effects. 

A much more complicated problem is the possible bias 
between DM and the galaxy distribution. The properties of 
large scale structures are moderately sensitive to the small 
scale bias (BBKS; Coles 1993; for review, Sahni & Coles 
1995), but available observations show that the spatial dis- 
tribution of DM and luminous matter can be biased on large 
scales as well. Indeed, while on one hand, in clusters of 
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Table 7. Parameters of RSE in mock catalogues in rcdshift space 

model runk b'^ Nthr ^rse frse Srse 

/i^^Mpc 

mocki 2.1 0.626 

mock2 2.0 0.586 

mocks 1.7 0.554 

mock4 1.9 0.563 



galaxies the observed ratio of galaxy and baryonic densi- 
ties is found to be Pgai/pgas ~ 0.2 (see, e.g.. White et al. 
1993), on the other hand, for example within Boots Void, 
pgai/pgas —> (Wcistrop et al. 1992). The existence of 'in- 
visible' structure elements, which are now seen as gas clouds 
responsible for weak Ly-Q absorption lines situated far from 
galaxies (~(5 — 6)/i^^Mpc, Morris et al. 1993; Stocke et al. 
1995; ShuU et al. 1996) can also be considered as an evidence 
in favor of a large scale bias. 

Such a large scale bias could be produced by the UV 
radiation from the first galaxy population during the re- 
heating of the universe (Dekel & Silk 1986; Dekel & Rees 
1987). Quantitative estimates (Demiaiiski & Doroshkevich 
1997, 1999a, b) show that it can increase the galaxy con- 
centration within the RSE by about a factor of 1.5 - 2. An 
indirect evidence in favor of such a bias was found in simu- 
lations by Sahni et al. (1994) as a suppression of structure 
formation in large regions around the maxima of gravita- 
tional potential, and by Doroshkevich, Fong & Makarova 
(1998) as a difference in the characteristics of spatial dis- 
tribution of the rare high peaks identified as 'galaxies' and 
of the main fraction of structure elements in simulated DM 
distributions. This bias can be essential for the reliability of 
discrimination between cosmological models. 

Unfortunately, such large scale bias cannot be simulated 
yet, since simulations of the galaxy formation in large boxes 
with the required resolution are impossible. This means that 
a bias can only be introduced by hand using simple plausible 
assumptions. Some such models as discussed by Cole et al. 
(1998) increase the 'galaxy' concentration within RSE. More 
detailed tests of these models will be discussed separately. 

Here we restrict our consideration to the analysis of 
simple mock catalogues prepared for the OCDM model. 
With the spatial and mass resolutions of our simulations 
(^IO^^Mq), we are compelled to identify 'galaxies' with 
selected DM particles. Four mock samples with different 
clustering properties were prepared and investigated. The 
simplest sample mocki was constructed with a threshold 
prescription depending on the local environmental density 
within a sphere of l/i~^Mpc around the particles. No 'galax- 
ies' are identified with particles with local density smaller 
than the mean density, and 'galaxy' tracers are selected ran- 
domly from the particles in overdense regions. The other cat- 
alogues use a smooth probability distribution proportional 
to the local environmental density within the same scale. 
The constant of proportionality was chosen to vary the de- 
gree of clustering of the mock samples, i.e. to get an autocor- 
relation function of simulated 'galaxies' in broad agreement 
with data (for the LCRS, cp. e.g. Tucker et al. 1997). The 
catalogues were normalized to the mean number density of 
galaxies, Ugai ~ 2 • 10~^h^ Mpc~^, that is equivalent to 
the observed galaxy density with the limiting magnitude of 
about Mr = -18. 
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Table 8. Proper sizes of RSE in mock catalogues 

{Nmem) (A) CTx {i^) 



model 
LCRS 
mocki-rd 
mock2-rd 
mocks-cm 
mocks -rd 
mock4-rd 



Tlnk 

2.1 
2.0 
1.9 
1.7 
1.9 



108 
119 
179 
174 
218 



25.3 
20.1 
19.3 
18.5 
18.3 
21.0 



6.0 

4.9 
4.4 
6.2 
4.5 
5.1 



12.3 
10.6 
10.5 

8.4 

9.4 

10.7 



3.1 

2.0 
2.1 
2.7 
2.2 
2.5 



5.7 
6.0 
5.9 
4.6 
5.1 
5.4 



1.3 

1.4 
1.4 
1.3 

1.4 
1.8 



\J-ypru)/ 

26.4 
23.0 
23.2 
22.8 
22.9 
27.5 



^ prvj 

1.4 
1.2 
0.9 
1.2 
0.9 
0.8 



\'-^ Sep) 

~38 
62 
53 
51 
49 
57 



-28 
57 
57 
53 
54 



All mean sizes and dispersions arc given by (6.2), (6.4) and Sec. 6.5 and are measured in h ^Mpc. 



The sample mocki with the threshold bias shows weakly 
enhanced clustering of the mock 'galaxies' with respect to 
the dark matter, the correlation function in redshift space 
shows a power law (r/ro)~^ with a correlation length ro = 
6.5/i~^Mpc and a slope 7 = 1.4. The correlation length of 
DM, also in the redshift space, is ro = 5/i^^ Mpc, and the 
slope is 7 = 1.3. 

The impact of the local environment allows to vary the 
clustering and to obtain mock samples with different prop- 
erties. Thus, the weakly clustered sample mock2 is similar to 
mocki (correlation length in redshift space ro — 6/i~^Mpc, 
slope 7 = 1.4) while mocks is intermediate (ro = 7/i~^Mpc, 
7 = 1.5), and mock4 strongly clustered (ro = 8/!.~^Mpc, 
7 = 1.6). For the sample mock2, we impose an additional 
threshold of no 'galaxy' identification for lower than mean 
density. We discuss such mock catalogues in our rectangular 
slices that allow a direct comparison with the DM catalogue 
studied above, ft demonstrates directly the influence of bias 
models. This is most important as a first step of the anal- 
ysis of the infiuence of bias. The selection criteria required 
for the more detailed comparison with observed catalogues 
as the LCRS (cp. the mock sample for ACDM presented in 
Fig. 1) will be imposed as a next step in a separate paper. 
The simple models described above reproduce only some fea- 
tures of the large scale bias and serve mainly as illustration 
of the potential of the used statistics for the bias discrimi- 
nation. More realistic models would be sensitive to the more 
broad density environment of particles, i.e. they are able to 
follow in more detail the expected interaction of large and 
small scale perturbations. Such models will be studied and 
discussed separately. 

The first step of our analysis repeats the approach uti- 
lized in Sec. 4 to define the sample of RSE with the required 
richness and overdensity. Such sets of RSE were found in 
redshift space for all four mock catalogues with parame- 
ters listed in Table 7. For all mock catalogues the thresh- 
old density of RSE was ~30% larger than that for the DM 
catalogue. For the mocki and mock2 catalogues, the small 
value of Nthr requires that the reproduction of a sufficient 
'galaxy' concentration within the RSE is also accompanied 
by a relatively low overdensity of RSE. Even so, the over- 
density is found to be at least two times larger than for the 
corresponding DM catalogue (see Table 2), and it is compa- 
rable with the overdensity found for the LCRS. For mocks 
and mock4 catalogues, the parameters of RSE are similar to 
each other, and to that found for the LCRS, and the over- 
densities reached are about three times larger than for the 
DM catalogue. For all mock catalogues the velocity disper- 
sions are the same as those found for DM RSE. These results 
alone show that the used models of bias provide an essential 
excess of 'galaxy' concentration within the RSE and, there- 



fore, they can be considered as a reasonable basis for the 
further more detailed investigation. 

The proper sizes of rich structure elements are found 
with the inertia tensor method, applied to the subsample of 
RSE. All proper sizes depend on the mass of the RSE. For 
all mock catalogues the scaling can be approximated by 



L = /.^-"^A, 



t = jy 



r U 

J^przu — ^ 



'D„ 



(8.1) 



0.3<U = Nraem/{Nraem) < 25. 



The mass-averaged length (A), width (ui), thickness (6) and 
the size of 'proto-wall', Dprw, are listed in Table 8 together 
with the corresponding dispersions. The stronger scaling 
found for the mock catalogues as compared to the DM cat- 
alogue is also caused by the biasing and therefore, this ap- 
proach can also be used to discriminate bias models. The 
main parameters of RSE listed in Table 8 are close to those 
listed in Table 4 for the OCDM model in redshift space. The 
shape and power indices found for the FDMST are also close 
to those listed in Table 5. As before, in comoving space the 
ID character of the 'galaxy' distribution dominates, i.e. the 
exponent pmst ~ 1, whereas in redshift space, the 'galaxy' 
distribution within RSE is similar to a 2D Poissonian dis- 
tribution with the exponent pmst ~ 2. 

These results show that the discussed methods reveal 
the influence of the used bias models. The impact of bias is 
also essential for the less massive structure elements, and es- 
pecially for the matter and 'galaxy' content in 'voids'. These 
regions will be further investigated with the MST technique 
in a separate study. 



9 SUMMARY AND DISCUSSION 

In this paper the properties of simulated spatial matter dis- 
tributions were studied for five cosmological models with 
CDM-like power spectra. The main parameters of simula- 
tions are listed in Table 1. The simulations were performed 
in large boxes in order to reproduce correctly the mutual 
interaction of large and small scale perturbations, and to 
obtain a representative sample of wall-like RSE. The broad 
set of considered cosmological models allows us to reveal 
the infiuence of main cosmological parameters on the forma- 
tion and evolution of the wall-like RSE, and to discriminate 
between these models. Our results show that the methods 
utilized in this paper are effective, and they yield a descrip- 
tion of the spatial matter distribution on large scales and, 
in particular, the characteristics of the RSE. 
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9.1 Main results 

The main results of our analysis can be summarized as fol- 
lows: 

(i) Simulations performed with the standard COBE nor- 
malized CDM-like power spectrum reproduce well both the 
wall-like RSE and the filamentary component of structures 
in LDR. Each component accumulates an essential fraction 
of matter, and it is equally important for the description of 
the joint network structure in the large scale matter distri- 
bution. 

(ii) The phenomenon of the strong matter concentration 
within the wall-like rich structure elements can also be repro- 
duced for suitable cosmological models. An essential fraction 
of DM, frse ~ 0.4, is compressed nonlinearly on the scales 
~(17— 25)/i~^Mpc that is less than the mean separation of 
these RSE by a factor of ~ 2 — 3. 

(iii) The RSE are usually disrupted into a system of high- 
density clumps that results in the growth of the inner veloc- 
ity dispersion. The rate and the degree of disruption depends 
on the parameters of the cosmological models. 

(iv) The comparison of observed and simulated parame- 
ters of the wall-like RSE allows us to discriminate between 
different cosmological models and to reveal the class of mod- 
els which can reproduce the main observed characteristics of 
the wall- like RSE. These are the ACDM model with flmh ^ 
0.15 - 0.25 and the OCDM model with O^/i ~ 0.25 - 0.35. 
Perhaps, promising results can be also reached for MDM 
models with similar parameters. 

(v) A large scale bias between the spatial distribution 
of DM and galaxies can increase the galaxy concentration 
within RSE by a factor of about 2 which essentially im- 
proves the simulated characteristics of RSE. The technique 
used above allows us to reveal reliably the influence of the 
biasing and to discriminate between bias models. 

(vi) The simulated parameters of DM structure are con- 
sistent to the theoretical expectations. The main cosmolog- 
ical parameters can be successfully reconstructed using the 
measured properties of the large scale matter distribution. 
After correcting for the bias and for selections effects, these 
methods can be applied to observed galaxy catalogues. 

(vii) Our results verify also the theoretical expectations 
with respect to the epoch of the RSE formation. At z = 1 
the fraction of matter accumulated by RSE with the chosen 
richness and overdensity drops by a factor ~2, and ai z — 3, 
it becomes negligible. 

The main statistical characteristics of the RSE are listed 
in Tables 2 - 5 in comparison with the properties found for 
the observed galaxy distribution. These results verify the ex- 
istence of a wall-like component with similar characteristics 
both in observations and simulations performed for suitable 
models. A similar range of cosmological models was recently 
separated by Cole et al. (1997) and by Bahcall & Fan (1998) 
from an comparison of observed and simulated properties of 
clusters of galaxies. The observations of supernovae at high 
redshifts (Perlmutter et al. 1998) are more consistent with 
the ACDM model with n„ ~ 0.3. 

Now there is some observational evidence of large mat- 
ter inhomogeneities at redshifts z ~ 0.5 — 1 and more 
(Williger at al. 1996; Cristiani et al. 1996; Quashnock et 
al. 1996, 1997; Connolly et al. 1996). Our analysis shows 



that for the considered models, these structures cannot be 
as common as at small redshifts. A more detailed statistical 
description of the absorption spectra of quasars is required 
to obtain the characteristics of structures at high redshifts. 
At small redshifts further progress in investigations of 
observed large scale matter distribution is linked with very 
large galaxy redshift catalogues as the 2dF redshift survey 
of 250,000 galaxies (CoUess 1998) and the million galaxy 
Sloan Digital Sky Survey (Loveday & Pier 1998). The avail- 
able surveys used above for the comparison with simulated 
structure parameters (the Durham/UKST Galaxy Redshift 
Survey, Ratcliffe et al. 1996; and the Las Campanas Red- 
shift Survey, Shectman et al. 1996) are not so representative 
and moreover they are limited in use to specific selection ef- 
fects (see discussion in Sec. 6.1 & 6.3). In spite of this, now 
and within the next few years results obtained with these 
surveys provide us with the best characteristics of observed 
large scale matter distribution. 

9.2 Methodical remarks 

The study of simulations is now the best way for the un- 
derstanding of the large scale matter distribution. The large 
boxes used for simulations allow us to obtain a representa- 
tive description of large scale perturbations and their mutual 
interactions with perturbations on smaller scales, as well as 
to obtain a representative statistic of RSE. Both factors are 
equally important for the successful reproduction of the mat- 
ter distribution observed in large galaxy surveys. 

The analysis of six simulations performed by Madsen 
et al. (1998) shows that simulations reproduce the theo- 
retical distributions only on the scales ~(0.1 — 0.15)Lbox- 
This means that realistic simulations of the observed large 
scale matter distribution is possible for Ltoa; > (400 — 
500)/i~^Mpc, whereas for smaller box sizes random varia- 
tions of parameters of large scale structure are expected. 

The analysis was performed mainly for rectangular 
slices with the size (500 x 500 x 50)/!."^Mpc^ which accu- 
mulate about 10% of the particles. Such an approach allows 
us to study a broad set of cosmological models with a rea- 
sonable precision and representativity. To test the possible 
impact of the selection used, the analysis was repeated for 
the full simulated sample of the ACDM model in comoving 
coordinates. The results are consistent with what was found 
above, and the difference is less than 10%. The comparison 
of results obtained for two ACDM models shows also the dif- 
ference between the structure parameters and the velocity 
dispersions ~10%. Variations of a^ei and as listed in Table 
1 are also ~10%. These results characterize the actual pre- 
cision reached in the investigation and shows that even for 
large boxes the main structure parameters are moderately 
sensitive to the random realization. 

The comparison of presented data with a similar sim- 
ulation (Cole et al. 1997) performed with higher resolution 
(P3M code) demonstrates the moderate dependence of the 
main simulated structure parameters on these factors. The 
properties of high density clumps are sensitive to the res- 
olution that distorts the FDMST for smaller lengths (in 
the comoving space). These distortions however disappear 
in redshift space. The same factor increases the simulated 
velocity dispersion. The main quantitative characteristics of 
the RSE are nonetheless sufficiently stable. 
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